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Abstract 

We introduce a new Large Eddy Simulation model in a channel, based on the 
projection on finite clement spaces as filtering operation in its variational form, for a 
given triangulation {7h}h>o- The eddy viscosity is expressed in terms of the friction 
velocity in the boundary layer due to the wall, and is of a standard sub grid-model 
form outside the boundary layer. The mixing length scale is locally equal to the 
grid size. The computational domain is the channel without the linear sub-layer of 
the boundary layer. The no slip boundary condition (BC) is replaced by a Navier 
(BC) at the computational wall. Considering the steady state case, we show that the 
variational finite element model we have introduced, has a solution {'Vh,Ph)h>o that 
converges to a solution of the steady state Navier-Stokes Equation with Navier BC. 

MCS Classification : 76D05, 76F65, 65M60, 

1 Introduction 

Numerical simulations of incompressible turbulent flows cannot be performed from the 
evolutionary Navier-Stokes Equations (NSE), 

, . 9tv + (v V)v - i/Av + Vp = f, 

^ V-v = 0, 

because of a great computational complexity due to the structure of the turbulence [30]. 
This is why various mathematical models derived from the NSE are used to simulate some 
features of turbulent flows, such as their statistical means or their large scales motions, 
this last way being known as "Large-Eddy Simulation" (LES), which is our concern in the 
present paper. 

LES has attracted much attention these last two decades, especially because of the increas- 
ing of computational ressources, enabling to enlarge the range of scales that LES models 
might simulate. Basically, LES aims at computing filtered fields such as v = G * v, G 
being a smooth transfer function [25, 30, 31, 32]. The filtering operation also might be 
carried out by solving PDE's [6, 7, 14, 16, 23, 24]. 
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Stresses that appear by filtering the non hnear term (v • V)v in the NSE, are considered 
to be diffusive, therefore often modeled by a turbulent diffusion term such as —V • (f^Vv), 
where ut is an eddy viscosity. One challenge of the LES is the determination of Vf 

In this paper, we study the case of a channel flow, periodic in the X1 — X2 axis for simplicity. 
The first idea is that the projection on finite element spaces, based on a given triangulation 
{Th}h>o-, is a natural filtering operation, so that we seek for v/j instead of v, where v/j is 
the projection of v on a suitable finite element space W/^. The second idea is that one 
can specifically model the eddy viscosity on the boundary layer by means of wall laws. 

Indeed, following Kolmogorov theory [ :], we consider the turbulence to be isotropic at 
scales small enough inside the flow domain. This assumption yields to take the eddy 
viscosity of a Kolmogorov-Prandtl-Smagorinsky form, vt = /i^|Z?v/j|^, h being the mixing 
length, which is the standard sub-grid model (SGM) [ ]. 

However, near the wall, turbulence is not isotropic and complexity is higher than far from 
the wall [29], so that standard SGM cannot be used there. Usual methods proceed as 
follows. 

First one uses the known structure of the turbulent boundary layer, as initially described 
by von Karman [ ] and fully developed by Schlichting [34]. Basically, the boundary layer 
may be split into two sub-layers, the linear sub-layer where the mean velocity profile is 
linear, and next, the log sub-layer where the mean velocity profile is specified by a log 
function. Notice that one can consider more sophisticated models to model the boundary 
layer [35], nevertheless always involving a log law. In all cases, those models involve an 
essential quantity which is the friction velocity u^, (see (2.4) in 2.1.ii below). 

Next, one splits the domain into two subdomains, the boundary layer, and the computa- 
tional domain which is the domain's part not containing the boundary layer. One then 
uses non linear boundary conditions at boundaries of the computational domain such as 
wall laws [26, 27, 28]. 

Based on the fact that today more computational resources are available to increase accu- 
racy for simulating the mean flow inside the log layer, we take as computational domain the 
domain's part without the linear sub -layer, using an eddy viscosity of the form vt = hu-^ 
inside de log layer, deduced from standard dimensional analysis [11, 26]. 

To conclude the modeling process, it remains to: i) specify how is calculated, ii) specify 
boundary conditions (BC) at computational domain boundaries, iii) fix the choice of the 
mixing length scale. 

i) We assume that log law holds inside the boundary layer. Thanks to invertibility of the 
non linear profil, we can define as m*(v,x), that satisfies suitable estimates (see 2. 2. iii 
and estimate (2.17)). 

ii) As the thickness of the linear sub-layer is very small compared to other scales involved 
in the problem, a Taylor expansion allows to deduce from the no slip condition at the flow 
domain boundary a Navier BC at the computational walls ( (3.6) in 3.1.i). This is as if 
the linear sub-layer would exert a friction over the log sub- layer. 

iii) The mesh yields natural numerical length scales hx , where hx is a diameter of any K G 
Tfi- Therefore, one takes of the form Vf = h'j^\DviT_\ on K G Th inside the computational 
domain, and ut = hKUi,{v,x) on G Th in the log layer (see 3.2.ii). 

1 DVh = (l/2)(Vv;, + Vv^) 
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Once this modeling process is completed, we get a model expressed in its variational form 
over finite element space W/j x M^, as described in 3.2.iv. So far as we know, this model 
is totally new, and can be generalized to more complex and realistic geometries thanks to 
a careful differential geometry analysis, which is a work under progress. 

We consider all over the paper the steady-state case, which is in coherence with the fact 
that in a permanent regime and for a developed turbulence, mean fields are steady, which 
is not in contradiction with the fact that fluctuations might be time dependent. 

We prove that this variational problem has a solution {vh,Ph) £ W/j x (Theorem 4.3) 
which converges to a solution (v,p) of the steady-state Navier Stokes Equation (NSE) with 
Navier BC (Theorem 4.4). 

The paper is organized as follows. We start with general setting. Then we derive from 
the NSE a description of the boundary layer, introducing the friction velocity. We specify 
the computational domain and Navier BC, and next we perform the finite element setting 
and get the model. Finally we state and prove Theorem 4.3 and Theorem 4.4. 
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work was performed. He also thanks E. Memin and the Fluminance team at INRIA 
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MTM2012-36124-C02-01. 



2 General framework 
2.1 Channel flow 

2.1.i. Geometry, equations and boundary conditions. Let be a channel periodic in the 
xi axis and X2 axis, of height 1 + 2(i in the xa-axis, for a small parameter d << 1, 

(2.1) S7/ = {x = (xi,X2,X3) G T2 X ]R3 s. t. -d<X3<l + d}, 

where T2 is the two dimensional torus defined by 

^ R2 , ^ 27r22 
T2 = — where T2 = — - — , 
h ^ 

and L > is a given length scale. Let Tf denote 



(2.2) 



/ 



{x G T2 X s. t. xs = -d 01 X3 = 1 + d}. 



The steady-state Navier-Stokes equations with the no-slip boundary condition are as fol- 
lows. 



(2.3) 



(v • V)v - uAv + Vp 
V- V 



f 





in Q 
in Qj ^ 
on Tf. 



The source term f is a body force per mass unit, typically the gravity. Assuming f G 
L^(r2/)^ = L^(il/), we know that this equation has a solution (v,p) G H^^'^/^($7)'^ x 
(see in [ ]), whose norms are bounded by constants that only depend on z^, 
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||f||o,2,0/ and d, andp is defined up to a constant. Uniqueness is known when |lf|lo,2,f7//i^^ 
is small enough. 

2.1.ii. Friction velocity. Let (v,p) be any solution of 2.3. We still denote by v the trace 
of V on Fj. We deduce from trace Theorems and Sobolev Theorem that v G W^'^{T f)^ = 
W^'^{T f). Therefore, it makes sense to consider Dv • n on Tj, where n denotes the 
outward-pointing unit normal vector at Tf, Dv = (l/2)(Vv + Vv*). We split the vector 
Dv ■ n into its tangential part and its normal part, 

(2.4) L>vn= (L>vn)^ + ((Dvn) -11)11. 
Let v^, G L^(Tf) be defined on Tj by 

(2.5) = t;.(v)(x) = {i^\{Dv • n),(x)|)^, 
called the friction velocity associated to v at x G Pj. 



2.2 Boundary layer description 

2.2.i. Length scale. Condition of uniqueness to system (2.3) is not satisfied in a steady- 
state turbulent regime. Let S be the set of solutions, which is a closed subset in L^(r2y). 
According to [ ], one can construct a probability measure /Li on 5. We consider the 
following velocity friction w^, G L^(rj) defined by, 

(2.6) w^ = j^v^iv)dfiiv). 
We finally define the meanfriction velocity by 

(2.7) = -^||it'*||o,2,ry- G II 

to which is associated the typical length scale A that characterises the boundary layer, 

(2.8) A = -, 
assuming 7^ 0. 

We conjecture that — 00 when !|f||o,2,ny^ 00. 

2.2.ii. Main assumption about the boundary layer structure. We focus on the bottom of 
J7j, {x3 = —d}, assuming that the boundary layer at the top {xs = 1 + d} has a similar 
structure. According to experiments (see in [-'> i]), we assume that in the boundary layer, 
the mean fluid velocity has a constant direction and only depends on the variable X3, 
which means v(x) = v{x3)e, for some fixed unit vector e. Without loss of generality, we 
can assume that e = ei is the unit vector pointing along the xi-axis. 

Notice that any plane P of the form P = {x^ = h} included in the boundary layer, and 
any vector N orthogonal at P being given, our assumption yields in particular v • N = 
at P. 

2.2.iii. Log law. Experiments and suitable assumptions about turbulence [11, 34] indicate 
that the boundary layer can be decomposed into two sub layers: 

• near the boundary where the velocity profile v is linear (linear sub layer). 
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• the next sub-layer specified by a log profile (log layer). 
To be more specific, we introduce the dimensionless variable 

(2.9) z+ = ^, 

and we consider the following continuous function defined on [0, by 

z+ if < z+ < 

where At ~ 0, 41 is the Von Karman constant. In practical calculations, one takes z^ ~ 20, 
and z^^r^ ~ 100, that measures the thickness of the logarithmic boundary layer, taken to 
be equal to lOOA. According to experiments [.i i], boundary layer thickness goes to zero as 
the Reynolds number goes to infinity. 

The profile v in the boundary layer at the bottom oiVtf is given by the formula 

(2.11) v{x?;) = u,l(^^ 

A similar description applies to the boundary layer at the top of $7/, {z = 1 + d}. 

2.2.iv. Friction velocity expressed as a function of the velocity. We still focus on the 
bottom. Any 2:3 > being given, Let 

(2.12) F{f3) = (3L{ap), a = ^. 

With this notation, equation (2.11) may be written as 

(2.13) v = F{u^), 
thanks to definition (2.8). 

Lemma 2.1. Let F : [0, +00) [0, +00) be defined by by (2.12). The function F is 
invertible, so that equation (2.13) can be written as u^, = F~^{v) at each given X3. 

Proof. We observe that the function L satisfies 

(2.14) hm ^ = Ci, 

x^0+ X 



(2.15) lim = C2, 

a;-i>oo log X 

where Ci and C2 are non-zero constants. As L is strictly increasing and continuous in 
(0, +00), then F is strictly increasing and continuous in (0,+oo). Also, by (2.14) F is 
continuous at /3 = with F{0) = 0. Moreover, by (2.15) lim F(x) = +00. Then F is 

bijective from [0, -l-oo) onto [0, +00), which yields the invertibility of F as claimed. □ 

Lemma 2.2. Denote h = F~^ . Then there exist a constant C = C(x3) > 0, bounded, 
such that 

(2.16) V7 > 0, /i(7) < Cixs) (1 + 7) < C(l + 7), C = sup C(X3). 
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Proof. Then h : [0, +00) H- [0, +00) is bijective and continuous. Also, 

lim SOLL = ii]XL — — = lim — - — - = lim — - = q. 

7-s>oo ^ t->oo F{t) t^-oo L{at) t^co log(at) L{at) 

The conclusion is a consequence of the continuity of /i. □ 

We deduce from Lemma 2.1, inequality (2.16) and because top and bottom layers have 
the same structure, that the friction velocity can be calculated at each x G BL from the 
velocity v, and satisfies the estimate 

(2.17) < = n^(v, x) < C{1 + |v|). 

3 Turbulence model 
3.1 Geometry and meshing 

3.1.1. Calculation domain. From now we assume that the boundary layer is included in 
the union of two strips, 

(3.1) BL = {-d < X3 < D/2 -d}U{l + d- D/2 < X3 < 1 + d}, 

where d < D « 1, d being the order of the linear sub layer, D the thikness of the global 
boundary layer. Standard numerical simulations are carried out in a sub domain of the 
flow domain that does not include the boundary layer at all, using a wall law [11, 27, 28] 
at artificial boundaries (walls). Our model includes the log layer, using a Navier BC based 
on a Taylor expansion as shown below. 

The computational domain is 

(3.2) = {x = {xi,X2,X3) G T2 X s. t. < X3 < 1} , 
the artificial wall being defined by 

(3.3) = {x G T2 X s. t. X3 = or X3 = 1}. 

3.1.ii. Boundary conditions. As in above, we focus on the bottom layer. By a Taylor 
expension we get 

dv 

(3.4) = v\x^=-d ^ v\x3=o - d-^\x3=o. 

From the view point of the domain ^l, v = VT-|r„, and d/dx-s = —d/dn at F^, where is 
the tangential part of v, defined by 

(3.5) V = + (v • n)n, 

by still denoting v the trace of v at F^, so far no risk of confusion occurs. Therefore, by 
remarks in 2.2.ii together with (3.4), we get 

d^T 1 

(3.6) v.n|r„=0, ^lr. = -^v., 
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which is a Navier boundary condition at the artificial wah, that expresses in some sense 
that the Hnear sub-layer exerts a friction on the log layer. Hence, system (2.3) becomes in 

n, 



(3.7) 



' (v • V)v - i/Av + Vp = f in n, 

V • V = in n, 

V • n = on r^o, 

-v^ on r,,,. 



dn d 

Navier-Stokes equations with Navier boundary conditions was studied before [1, 2, 3, 10, 
37], and existence of a solution to system (3.7) is already ensured. 

3.1.iii. Variational formulation. Let us define the spaces 

W(l^) = {w G Hi(l^), wn|r„=0}, 



M{Vl) = {qeL^{Vl), [ qdx 

Jnc 



0}, 



by reminding that H^(il) = i7^(il)^. Strictly speaking. The space M{Q,) is isomorphic to 
the quotient space L'^{fl)/R, endowed with the usual quotient norm 

(3.8) ||p||M = inf |b||o,2,n. 

It also may be viewed as a closed subspace of L'^{Q) endowed with the L^(r2) norm. 

The space W(J7) is endowed with the norm, denoted || • ||i,2,Q- As a consequence of 
Korn's inequality, the following usefull estimate holds, 

(3.9) VvGW(l^), ||v||i,2,n <C(||L'v||o,2,n + ||v||o,2,r„), 
of proof of which being carried out in [10]. 

Let a, b and G the forms defined by 

(3.10) a(v,w) = u{Dv,Dw)n, 

(3.11) Kz;v,w) = i[((z-V)v,w)o-((z-V)w,v)n] 

(3.12) G(v,w) = ^(v„wOr„, 

for z, V, w € H^(f2). Recall that when z, v, w G W(J7) and V • z = 0, then 6(z; v, w) = 
((z • V) V, w)q, and (Vz, Vw)q = {Dz, Dw)^. Also remark that when v G b, then v = Vt- 
at . 

We say that a pair {v,p) G W(r2) x Af (il) is a weak solution of the boundary value (3.7) 
if it satisfies 



(3.13) 



b{v;v,w) + a{v,w) - {p,V ■ w)q + G{v,w) = 

{V■^r,q)n = 0, 



for any (w, g) G W{n) x M(17). 

3.1.iv. A priori estimate and existence result. Assume f G W(i7)'. Let (v,p) be any solu- 
tion of Problem (3.13), and take v = w in (3.13). From the standard formula 6(v; v, v) = 
that holds since V • v = and v • n = at F^^,, we get 

(3-14) i^||-Dv||o,2,Q + ^||v||o,2,r„ = (f, v) , 
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from where we deduce 



(3.15) ||v[|i,2,n < ^||f||w{n)', k = mm , 
by using (3.9). 

3.2 Finite element setting 

3.2.i. Triangulation. Let D C denotes the sample box D = [0, L]^ x [0,1]. The 
computational domain Cl may be viewed as the periodic reproduction of D in the xi — 
X2 axes. Let {Th}(h>o) be a regular familly of triangulation of D, compatible with the 
periodicity of the domain: The restriction of the grid to the planes xi = and xi = D 
is the same, and the restriction of the grid to the planes X2 = and X2 = D is the same. 
Reproducing this triangulation by periodicity, we get a regular triangulation of fi, still 
denoted by {Th}(h>o)- 

In the following, for each K € Th, hx = diam{K) denotes the diameter of A', and 

(3.16) h = max hx- 

3.2.ii. Eddy viscosities. We assume isotropy of the turbulence inside the domain defined 
by Qin = 0,\BL, the boundary layer BL being defined by (3.1). This yields to consider in 
Qin the eddy viscosity ut^in to be of the same form as in usual Sub-Grid-Models of Prandtl- 
Kolmogorov-Smagorinsky type, where following [Ll], we take in each K the length scale 
equal to hx, leading to consider i't,in to be of the form 

(3.17) ut,M = Cl li^l^vl, 

Cs > being an empirical constant, Ia denotes the characteristic function for any set A. 

In the boundary layer part, = BL fi Vt, turbulence is no longer isotropic and depends 
on the friction velocity. Taking again hx as typical length scale and by a dimensional 
analysis argument [11], we define the eddy viscosity vt,w in by 

(3.18) t't,«,(v)=C^ ^ /liC lii"u*(v,x), 

where Cw > is an empirical constant and Uj, is expressed in 2.2.iv. 
Finally, the eddy viscosity we consider is of the form 

(3-19) vt = i^tM = lf7,„ft,m(v) + ln^i^t,w{'y), 

3.2.iii. Finite element spaces. The model is a mixed formulations, based upon pairs of 
finite element spaces (W/i,M/i) C W(i7) x M(r2), associated to the family of regular 
triangulations {Th}h>o of O in the sense of Ciarlet [I;']. We assume that the family of 
pairs of spaces {(Wh, Mfi)}h>o satisfies the following hypothesis: 

Hypothesis 1: The family of spaces {Wh x Mh}h>o is an internal approximation of 
W(Q) X M(Q): For ah (w,p) G W(r2) x M{^}) there exists a sequence {{vh,Ph)}h>o such 
that {■Vh,Ph) e X Mh, and 

lim(||v - Vh\\i,2,n + \\p - Ph\\o,2,n) = 0. 

h—>-0 
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Hypothesis 2: The family of pairs of spaces {(W/j, M/i)}/i>o satisfies the uniform discrete 
inf-sup condition : There exists a constant a > such that 

rn nn\ II II ^ (V-W/i,g/i)Q 

(3.20) a\\qh\\o,2,n< sup — ^ ^ , for all g,, e M,, 

There is a wide literature about finite element spaces satisfying those properties (Cf. [•'")], 

[9], im 

3.2.iv. The model. Our LES model is expressed by the following variational problem: 
Find {vh,Ph) G W/, x Mh such that for all iwh,qh) G x Mh, 



(3.21) 



^(v/i; v/i, w/,) + a(v,i, w?,) + c(v?,; w/i) + 

G(v/i, w/i) - (p/i, V • w/i)n = (f,w/i), 
(V-Vft,g/i)n = 0; 



the form c being defined by 

(-3 22) '^(^; = (^*(^) or 

c(v;w) = (i/t,j„(v)L>v,Dw)n,„ + (t't,«,(v)93V,a3w)n„. 

The second expression neglects the tangential eddy viscosity in the boundary layer, which 
is very small compared to the normal one. 

4 Analysis of the model 
4.1 Technical results 

We state in this sub-section some technical results concerning the eddy viscosities and the 
associated turbulent diffusion form c, that are needed by our analysis. 

4.1.i. eddy viscosties estimates. 

Lemma 4.1. There exists a constant C > depending only on the aspect ratio of the 
family of triangulations {Th\h>o such that 

(4.1) pt{'Vh)h,oo,n<Ch^'^\\ v/i||i,2,fi5 for all v/j € W/j. 

Proof. We start with the internal part of the eddy viscosity ut,in- Consider v/^ G W/^. As 
Vv/i is piecewise continuous, there exists K G Th such that 

\Wt,ini'yh)\\o,oo,n = \Wt,ini'yh)\\o,oo,K < Cs^K II Vv/i ||o,oo,_ft:. 

By a standard finite element inverse estimate (Cf. [4]), 

—3/2 1 1 

<CV llVvhIl 

for some constant C > depending only on the aspect ratio of the family of triangulations. 
Then, 



(4.2) \h,ini^h)h,oo,n < CClh^-^'^ ||Vv;,||o,2,ic < CClh^l^ II Vv 



h\\o,2,n, 
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Next, we analyze the wall eddy diffusion i^t.w There exists some element K G Th such 
that 

lkt,«)(v/i)||o,oo,n = \\l't,w{'^h)\\o,oo,K < Cw hx (1 + IIv/jIIo.oo,^), 

where in the last inequality we have used (2.5). Using the inverse estimate (Cf. [ ]), 

— 1/2 

||v/i||o,oo,x <Chj^' ||Vv/i||o,2,x we deduce 

K)|| 0,2, for some constant C > 0. 

Combining this estimate with (4.2) and ||Vv/j||o,2,r2 ^ ||v/i||i,2,0) (4.1) follows. □ 
4.1.ii. Turbulent diffusion operator properties. 

Lemma 4.2. The form c defined by (3.22) satisfies the following properties: 
i) c is non-negative, in the sense that 

c(v;v) > 0, for all v G H'^{nf. 

a) Assume that the family of triangulations {Th\h>o is regular. Then, for any Vh, w/j G 
(4-3) |c(v,,;w,,)| < C/i^/^ ||v/i||i_2,nl|wft||i,2,o, 

for some constant C > depending only on d, and the aspect ratio of the family 
of triangulations. 

Hi) Assume that the family of triangulations {Th\h>o is regular. Let {^h}h>o o-nd {'Wh}h>o 
be two sequences such that v/j, G W/j. Then, if both sequences are bounded in 
Hi(0)^ 

(4.4) limc(v,,;wft) = 

Proof. 

i) Let V G Hi(J7). Then, 

c(v;v)= / i^t(v) iDvpdx > 0. 
Jn 

ii) By estimate (4.1), 

\c{vh;Wh)\ < ||ft(v/x)||o,oo,f7 ||v/i||l,2,n||w/i||i,2,r2 
^ 1.1/2 11 |i2 II II 

< on' \\^h\\i^2,n\\^h\\i,2,n- 

iii) Statement (4.4) directly follows from (4.3). 

□ 
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4.2 Existence result 

Problem (3.21) is a set of non-linear equations in finite dimension. These non-linearities 
are due to several effects: the convection operator, the eddy viscosity, and the wall-law 
boundary conditions. The space W(ri) is a closed sub-space of H^(Q). Our main result 
is the following. 

Theorem 4.3. Let {Th}h>o be a regular family of triangulations of the domain il. Let 
{(W/i, M/i)}/i>o he a family of pairs of finite element spaces satisfying Hypotheses 1 and 
2. Then for any f G W(r2)' the variational problem (3.21) admits at least a solution, that 
satisfies the estimates 

(4.5) ||v/,||i,2,Q < 

(4.6) \\Ph\\o,2,n < 

where C > is a constant depending only on d, 17 and the aspect ratio of the family of 
triangulations. 

Proof.. We prove the existence of solution in two steps. 

4.2.i. Step 1: Existence of the velocity. Let us define the mapping : — > as 
follows: Given Zh £ W/i, 

for any Wh G W/j. This equation has a unique solution as its r.h.s. defines a linear bounded 
functional on W/j. Moreover, the functional is continuous as all functions that appear 
in its definition are continuous on the finite-dimensional space W/j. 
Consider the sub-space Zh of defined by 

Zh = {w/i G Wh such that (V • Wh, Qh) = 0, for all qt G Mh }. 

Zh is a non-empty closed sub-space of H^(r2). Then it is a Hilbert space endowed with 
the H^(r2) norm. Let G Z^. Then, as 6(z/i; z/^, z/^) = and c is non-negative, 

(^/i(z/i),Zh) > a{'^h,'^h) + G{zh,Zh) - {{,Zh) 

> V ||-D(zfe)||o,2,n + ^l|zh||o,2,r™ - l|f|lw(Q)'l|z/i||i,2,n 

Ck 2 ll^llw(n)' 

> -Yl|zh||i,2,n- 2C7k 

where we have used (3.9) and Young's inequahty. We deduce 

||f 

(4.7) Vzft, G Zh such that ||z/,||i,2,f7 = — , {^h{^h).'^h) miO) > 0, 

Consequently, by a classical variant of Brouwer's Fixed Point Theorem (Cf. [ ]), the 
equation 

(4.8) b{vh;\rh,y^h) + a{\rh,\Vh) + c{^rh;wh) + G{vh,v^h) = {i,'^h) ^v^h^Zh 

ll^ilwfn)' 

admits a solution Vh G Zh such that ||vft||i^2,n < — , which precisely is (4.5) by 

changing C in C^^. 



Ck ""^llf I|w(q)') K = min(z^, 



d 



C7K-i||f||w(fi)' U-'m\w(ny[l + h^'^] + + ^ + l) 
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4.2.ii. Step 2: Existence of the pressure. Let the operator ■ Mf^ i— t- defined by 
\/qh G Mh, {Qh{qh),^yi) = (V • ^rh,qh)n, for all e W/^. 

Then = Im{Qh)^. As Im{Qh) is closed, then = Im{Qh)- As v/j is a solution of 

(4.8) , then $/j(v/i) G Z^. Consequently, there exists some discrete pressure ph such that 
(^/i(v/i), w/j) = (V • Vh,Ph)n, for ah w/j € Wh- Thus, the pair {^h,Ph) solves problem 
(3.21). The estimate for the norm of the pressure is obtained via the discrete inf-sup 
condition (3.20), 

for some constant a > 0. By estimates (4.3) and some standard estimates, 

|v/»||?,2,n(l + Ch^/^) + U\\^^h\\l,2,n (^1 + ^) 
+ l|f|lw(n)' l|w/i||l,2,n- 

Then, the pressure estimate (4.6) follows from the velocity estimate (4.5). □ 
4.3 Convergence 

We now prove the convergence of the solution provided by method (3.21) to a weak solution 
of the Navier-Stokes boundary value problem model (2.3). 

Theorem 4.4. Under the hypotheses of Theorem 4-3, the sequence {i'Vh,Ph)}h>o contains 
a sub-sequence strongly convergent in H^(r2)^ x L^(r2) to a weak solution (v,p) S W(ri) x 
Lq(J7) of the steady Navier-Stokes equation (2.3). If this solution is unique, then the whole 
sequence converges to it. 

Proof. The proof is divided into 7 steps. 

4.3.i. Extracting sub sequences. By estimates (4.5) and (4.6), the sequence {{vh,Ph)}h>o 
is bounded in the space W(f2) x Lq(Q) which is is a Hilbert space. Therefore, this sequence 
contains a subsequence, that we denote in the same way, weakly convergent in W(r2) x 
Lq{Q) to some pair (v,p). As the injection of H^(Q,) in L'^(il) is compact for 1 < < 6, 
we may assume that the subsequence is strongly convergent in L''(il) for 1 < q < 6, and 
so in particular in L^(il). 

Also, the injection of i^^^'^{Ty^) into L^(r^) is compact. Then we may assume that the 
sequence {v/i, }h>o is strongly convergent to vi in L^(r^). 

4.3.ii. Taking the hmit in the diffusion terms. Let (w, q) € W(0) x Lq{Q). By Hypothesis 
1, there exists a sequence {{'Whj Qh)}h>o such that {wh,Qh) £ x which is strongly 
convergent in H-'^(O) x L'^{Q) to (w, g). 

As a is bilinear and continuous, 

(4.9) lim a{vh,Wh) = a(v, w). 

Next, since the sequences {v/i}/j>o and {'Wh}h>o are bounded in H-^(ri), we deduce from 
Lemma 4.2, 

(4.10) limc(v/,;w/,) = 0. 
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Finally it is straightforward to check that 

(4.11) limGK;wO = 0. 

4.3.iii. Taking the limit in limit in the convective term. We have 

|(vh • Vv/,, w/,)n - (v • Vv, w)n| < 
|((vh - v) • V\rh,Wh)n\ + |(v • V(v/j - v),w)n| + Kv • Vvh,Wh - w)nl < 

||vh - v|lo,4,Q ||Vv/i||o,2,r2 |jw/i||o,4,Q 

3 

+ ^ \{dj{vhi - Vi),VjWi)n\ + \\v\\o,i,n ||Vv/,||o,2,n llw/^ - w||o,4,n, 
«,i=i 

where we denote Vh = (w/ii, 'y/i2) I'/is)- All terms in the r.h.s. of the last inequality vanish in 
the limit because {v/j}/i>o is strongly convergent in L^(ri), {diVhi}h>o is weakly convergent 
in L-^(J7) and {w/j}/i>o is strongly convergent in H^(il). Then, 

(4.12) lim((v;, • VvO,wOc = ((vVv),w)n. 

/i-i-O 

Similarly, lim((v/j • V)wh,Vh)n = ((v • Vw),v)n, and then 

lim b{\'h; ^rh,Wh) = &(v; v, w). 

4.3.iv. Taking the limit in the pressure terms. Since {V • v/i}/i>o is weakly convergent in 
L^(r2) to V • and {qh}h>o is strongly convergent in L^(Q) to q, 



Finally, we obviously have 



lim (V • V,,, qh)n = (V • v, q)n. 



lim {ph, V • Wh)n = (p, V • w)o. 

/i->0 



Consequently, the pair (v, g) is a weak solution of Navier-Stokes equations (3.13). 
4.3.V. Strong convergence of the velocities. Set = in (3.21). Then 

\\Dvh\\l,2,n + '^hh\\o,2,r^ = (f,v/,) -c(vft;v/,). 

By Lemma 4.2 iii), lim c{vh',Vh) = 0. Therefore, 

/i->0 

£o 11-^'^'^ 1 1 0.2,0 + ^l|vh||o,2,r„) = (f, v) = 1/ lll^vllg^a.n + ^l|v||o,2,r„> 

where the last equality occurs because (v, q) is a weak solution of Navier-Stokes equations 
(3.13). As W(r2) is a Hilbert space and {'Vh}h>o is weakly convergent to v, this proves 
the strong convergence, since 



(4.13) (^u \\Dw\\l2,n + ^ l|w||o,2,r») 

is a norm equivalent to the H^(r2) norm by (3.9). 
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4.3.vi. Strong convergence of the pressures. We use the discrete inf-sup condition to 
estimate \\ph — p\\o,2,Q- There exists a sequence {Ph}h>o such that Ph £ for all /i > 
which is strongly convergent in Lq(Q) to p. We shall show that lim \\ph — Ph\\o2n = 0. 

Let w/i G W/j. We have 

{Ph-Ph,^-^h) = Kv/i; v/i, w;,) - 6(v; v,w?i) + a(v;, - v,Wft) +c(v/i; w/j) 
+ G{wh-^,^h) + (p- Ph,^ -^h)- 

As 

Kv/i; V/,, w/,) - &(v; v,vi^ft,) = 6(vft,;v/, - v,w/,) + 6(v/i - v;v,w/i) 

< C ||v/j - v||i,2,n (||v/i||l,2,Q + ||v||i,2,f7), 

using (4.3) and the continuity of a we deduce 

{Ph - Ph,^ ■ W/i) < C [\\wh - v||i,2,Q (||vh||i,2,n + ||v||l,2,n) + T^\\D{^/h - v)||o,2,n 
V 



+ ^^^^ l|v/i||?,2,n + ^ - v||o,2,r„ + Hp - -P/i||o,2,n 



Wft||i_2,r2. 



As lim ||v/j - v||i 2 n = 0, then by Hypothesis 2, lim \ph - -Ph||o,2,Q = 0. Then ph strongly 
converges to p in 

4.3.vii. Uniqueness. It remains to prove that if the Navier-Stokes equations (2.3) admit 
a unique solution (v,p), then the whole sequence {(v/j,p)}/i>oconverges to it. This is 
a standard result that holds when compactness arguments are used, which is proved by 
reductio ad absurdum: Assume that the whole sequence does not converge to {v,ph). Then 
there exists a sub-sequence of {(yh,Ph)}h>o that lies outside some ball of W(r2) x Lg(ri) 
with center {v,p). Then the preceding compactness argument proves that a sub-sequence 
of this sub-sequence would converge to the unique solution (v,p), what is absurd. □ 



References 

[1] Beirao Da Veiga, H.: On the regularity of flows with Ladyzhenskaya shear-dependent 
viscosity and slip or nonslip boundary conditions, Comm. Pure Appl. Math. 58 
552577 (2005). 

[2] Beirao Da Veiga, H.: Vorticity and regularity for flows under the Navier boundary 
condition. Commun. Pure Appl. Anal. 5 907918, (2006) 

[3] Berselli, L.C.: An elementary approach to the 3D Navier-Stokes equations with 
Navier boundary conditions: existence and uniqueness of various classes of solu- 
tions in the flat boundary case. Discrete Contin. Dyn. Syst. Ser. S, 3 (2), 199-219, 
(2010) 

[4] Bernardi, C, Maday, Y., Rapetti, F.: Discretisations variationnelles de problemes 
aux limites elliptiques. Mathematiques &; Applications, 45. Springer- Verlag, Berlin 
(2004) 

[5] Bernardi C. and Raugel, G.: Analysis of some finite elements for the Stokes problem. 
Mathematics of Computation, 44 (169), 71-79 (1985) 



14 



[6] Berselli, L., Lewandowski, R.: Convergence of approximate deconvolution models to 
the mean Navier-Stokes Equations. Annales de I'lnstitut Henri Poincare (C), Non 
Linear Analysis, 29, 171-198, (2012) 

[7] Berselli, L. C, Iliescu, T., Layton, W. J.: Mathematics of Large Eddy Simulation of 
Turbulent Flows. Springer, Berlin (2006) 

[8] Brezis, H.: Analyse fonctionelle: Theorie et applications. Masson, Paris (1983) 

[9] Brezzi, F., Fortin, M.: Mixed and Hybrid Finite Element Methods. Springer- Ver lag, 
New York (1991) 

[10] Bulicek, M., Malek, J., Rajagopal, K. R,. Navier's slip and evolutionary Navier- 
Stokes-like systems with pressure and shear-rate dependent viscosity. Indiana Univ. 
Math. J., 56 (1), 51-85, (2007) 

[11] Chacon ReboUo, T., Lewandowski, R.: Mathematical and numerical foundations of 
turbulence models and applications. To appear in Springer, (2013) 

[12] Ciarlet, Ph.: The Finite Element Method for Elliptic Problems. Classics in Applied 
Mathematics, Vol. 40. SIAM, Philadelphia (2002) 

[13] Clement, P.: Approximation by finite element functions using local regularization. 
RAIRO Anal. Numer. 9 (R-2) 77-84 (1975) 

[14] Germane, M.: Differential filters for the large eddy numerical simulation of turbulent 
flows. Phys. Fluids 29 (6), 1755-1757 (1986) 

[15] Girault, V.; Raviart, P. A. : Finite element methods for Navier-Stokes equations: 
theory and algorithms. Springer- Verlag (1986) 

[16] Foias, C, Holm, D.D., Titi, E.S.: The Navier-Stokes-alpha model of fluid turbulence. 
Physica D, 152, 505-519, (2001) 

[17] Hughes, T. J. R., Feijoo, G., Mazzei, L:, Quincy, J-B: The variational multiscale 
methoda paradigm for computational mechanics. Comp. Meth. Appl. Mech. Enrgrg. 
166 (1-2), pp. 3-24 (1998). 

[18] Hughes, T. J. R., Mazzei, L., Jansen, K. E.: Large eddy simulation and the variational 
multiscale method. Comput. Vis. Sci. 3 47-59 (2000) 

[19] Hughes, T. J. R., Mazzei, L., Oberai, A., Wray, A.: The multiscale formulation of 
large eddy simulation. Phys. Fluids 13 (2), 505-512 (2001) 

[20] Hughes, T. J. R., Oberai, Mazzei, L.: Large eddy simulation of turbulent channel 
flows by the variational multiscale method. Phys. Fluids 13 (6), 1784-1799 (2001) 

[21] von Karman, T. Mechanische Ahnlichkeit und Turbulenz. Nachr. Ges. Wiss. 
Gottingen, Math. Phys. Klasse, 58, (1930). 

[22] Kolmogorov, A.N.: The local structure of turbulence in incompressible viscous fluids 
for very large Reynolds number. Dokl. Akad. Nauk SSR, 30, 9-13, (1941) 

[23] Layton, W., Lewandowski, R.: Analysis of an eddy viscosity model for large eddy 
simulation of turbulent flows. Journal of Mathematical Fluid Dynamics, 4, 374-399, 
(2002) 



15 



[24] Layton, W., Lewandowski, R.: A simple and stable scale similarity model for large 
eddy simulation: energy balance and existence of weak solutions. Applied Maths. 
Letters, 16, 1205-1209, (2003) 

[25] Lesieur, M., Metais, O., Comte, P. Large-eddy simulations of turbulence. Cambridge 
University Press, (2005) 

[26] Lewandowski, R.: Analyse Mathematique et Oceanographie. Masson (1997) 

[27] Mohammadi, M., Pironneau, O.: Analysis of the fe-epsilon turbulence model. Masson, 
Paris, (1994) 

[28] Pares, C: Existence, uniqueness and regularity of solution of the equations of a 
turbulence model for incompressible fluids. AppHcable Analysis 43, 245-296 (1992) 

[29] Piomelli, U., Balaras, E.: Wall-layer models for large-eddy simulations. Annu. Fluid 
Mech. 34, 349-374, (2002) 

[30] Pope, S.B.: Turbulent flows. Cambridge University Press, (2000) 

[31] Rodi, W., Fcrzigcr, J.H., Brcuer, M., Pourquic, M.: Status of largc-cddy simulation. 
Workshop on LES of flows past bluff bodies (Rottach-Egern, Tegernsee, 1995). ASME 
J. Fluid Engng., 119, 248-262, (1997) 

[32] Sagaut, P.: Large Eddy Simulation for Incompressible Flows. Springer, Belin (2002) 

[33] Scott, R., Zhang, S.: Finite element interpolation of non-smooth functions satisfying 
boundary conditions. Math. Comput. 54, 483-493 (1990) 

[34] Schlichting, H. Boundary Layer Theory, 8th. Edition. Springer (2000) 

[35] Spalding, D.B.: A single formula for the law of the wall. Transactions of the ASME, 
Series E: Journal of Applied Mechanics, 28, 455-458 (1961) 

[36] Temam, R.: Navier-Stokes equations. AMS Chelsea (2001) 

[37] Verfiirth, R.: Finite Element approximation of incompressible Navier-Stokes equa- 
tions with slip boundary conditions. Numer. Math. 50, 697-721 (1987) 



16 



